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Abstract 

In this paper, we develop an analytical approach to predict the field-induced 
alignment of cantilevered graphene nanoribbons. This approach is vali- 
dated through molecular simulations using a constitutive atomic electrostatic 
model. Our results reveal that graphene's field-oriented bending angle is 
roughly proportional to the square of field strength or the graphene length 
for small deformations, while is roughly independent of graphene width. The 
effective bending stiffness and the longitudinal polarizability are also found 
to be approximately proportional to the square of graphene length. Com- 
pared with carbon nanotubes, graphene nanoribbons are found to be more 
mechanically sensitive to an external electric field. 

1. Introduction 

Graphene's electronic gap tunable in external electromagnetic fields [U 
El [31 H] makes it promising for a number of potential applications in na- 
noelectronic devices [HI El [7J. Since graphene is often supposed to work in 
a transverse electric field in such devices and its electronic transport prop- 
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6 erties strongly depend on its atomic structure [SJ HJ [TO], understanding of 

7 graphene's mechanical behaviors in an electric field is of great importance for 
s nanoelectromechanical systems based on graphene. However, as a novel is- 

9 sue, the mechanical response of graphene to an external field has not yet been 

10 reported up to date. How does graphene deform in response to applied elec- 

11 trie fields? To answer this question, we developed a simple model to predict 

12 the field-induced alignment of cantilevered graphene nanoribbons (GNRs), 

13 demonstrating the coupling between the graphene's field-induced bending, 

14 molecular stiffness and electric polarization. This model is validated through 

15 molecular simulations, in view of the difficulties of experimental quantifica- 

16 tion of this electromechanical effects in nanoscale. 

17 If a thin nanostructure is brought into an electric field, electric polariza- 

18 tion effects will induce a moment of force, which tends to orient the nanos- 

19 tructure toward the field direction [TT]. This moment will make the nanos- 

20 tructure bent if one end of the nanostructure is fixed on a substrate (see 

21 Fig. [T|. This alignment has been observed on carbon nanotubes (CNTs) in 

22 early experiments [12], and then exploited in designing nanorelays [13] and 

23 field emission devices [2] . As can be expected in view of the large similarity 

24 in their atomic structures and mechanical properties [15], graphene nanorib- 

25 bons (GNRs) should exhibit similar mechanical behaviors in an electric field. 

26 Furthermore, unlike CNTs, graphene exhibits strong mechanical anisotropy 

27 in its transverse direction, with a lateral stiffness about 30 times lower than 

28 that in the longitudinal direction [161 HZ]- This high lateral structural flexi- 

29 bility makes GNRs ideal field-sensing materials in resonators [5], transistors 

30 [T5] or sensors [T5] . 
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Figure 1: Alignment of a cantilevered GNR (L rs 20nm, periodic in width direction) to an 
electric field E = 0.3V/nm. The color scale of atoms represents the intensity of induced 
charges q. 

31 2. Simulations 

32 In this work, molecular simulations were performed to compute equilib- 

33 rium structures of cantilevered GNRs in an electric field, using an energy 

34 optimization method (full details of the simulation techniques can be found 

35 elsewhere [201 EI])- The principle of this method is to minimize the potential 

36 energy of each atom, which consists of an internal potential due to the C-C 

37 chemical bonds and the long-range interaction, and of an external potential 

38 induced by applied electrostatic fields. The internal potential is calculated 

39 using the adaptive interatomic reactive empirical bond order (AIREBO) po- 
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40 tential function [22], which has been used in recent studies on mechanical 

41 properties of CNTs [23] and GNRs [24] . The external potential describes 

42 the electrostatic interaction between the charges, the dipoles and the ex- 

43 ternal field, it is computed using an atomic charge-dipole interacting model 

44 [251 127] 128] . which has been validated through charge- injection experiments 

45 using an atomic force microscope [29|. 

46 3. Modeling 

47 The main reason of the field-induced alignment of a GNR is the effect 

48 of electric polarization, by which a quantity of positive and negative charges 

49 are shifted to opposite directions in graphene (see Fig. [TJ. A moment of 

50 a force pair arises from the electrostatic interactions between the field and 

51 the induced charges, and makes the graphene bent into the field direction. 

52 This field-driving moment highly correlated with the GNR's polarizability is 

53 resisted by the mechanical lateral stiffness of graphene, due to the repulsive 

54 interactions between 7r electrons and the rotation of a bonds. Considering 

55 the correlations between the electric polarizability, the bending stiffness, and 

56 the geometry of GNRs, the curvature uj of a cantilevered GNR in an electric 

57 field can be calculated as follows: 

_ M _ E 2 {a//-a ± ) sin [2(0 O - 0)] 

58 where M is the bending moment induced by the electric polarization J3U] , 

59 K denotes the effective bending stiffness of GNR, L stands for the graphene 

60 length, 9 is the deflection angle of GNR, E = \E\, and 9q represent the 
ei strength and the direction of the electric field, respectively, a/j and a± stand 
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62 for the longitudinal and the transverse molecular polarizabilities of GNR, 

63 respectively. 9 is defined as the angle between the initial axis of graphene 

64 and the vector from one graphene end to another after deformation. Since 

65 the curvature u can be approximated as to = 29 /L and usually a// » a±, 
ee the governing equation can therefore be written as 



LE 2 a*u sin [2(0 O 



(2) 
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Figure 2: K* vs L from simulations, in which the periodic condition was applied in width 
direction. Symbols show simulation data and the dashed line stands for a fitting curve. 
The inset shows the main stress distribution in a curved graphene. 



67 where a*^ and K* stand for the quantities of a// and K per unit width, 
es respectively. In Eq. |2j the parameter representing the mechanical resistance 
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of a GNR is the effective bending stiffness K. Here its values were directly 
computed from atomic simulations by applying a mechanical force at the 
free end of the GNRs in absence of electric field , instead of using a conven- 
tional formula for macroscopic continuous media as the product of the elastic 
modulus and the area moment of inertia, in order to avoid the problem of 
definition of the wall-thickness of a one-single-atom thick layer |31j . Note 
that previous studies showed that the bending stiffness of a CNT is an inde- 
pendent parameter not necessarily related to the representative thickness by 
the classic formula [32J . Our simulation results show that K* is about 20eV 
when L < 20nm and is roughly proportional to L 2 when GNRs get longer 
(see Fig. [2]). These simulation data can be fitted using a simple equation as 

K* = AL 2 (L > 20nm) (3) 

where A = 0.35eV-nm -2 for GNRs with either armchair or zigzag edges. 
For comparison, we also calculated the effective bending stiffness of CNTs. 
It is found that the value of K of a GNR (L pa lOnm, K pa 19.2eV-nm) 
is about 20 times smaller than that of a (5, 5) CNTs of the same length 
(K pa 570eV-nm). This large difference in the transverse (lateral) stiffness 
implies that the alignment of GNRs can be much more significant than that 
of CNTs for a given magnitude of electric polarization. 

Another important parameter in Eq. |2jis a*/. Its value was determined 
using electrostatic simulations based on the atomic charge-dipole model [25J . 
In these simulations, the value of ol* j was calculated from the definition 
p = Ea, where at is the the matrix of molecular polarizability, p and E 
stand for the vectors of the induced molecular dipole (see Fig. [3] (a)) and the 
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(b) L(nm) 

Figure 3: (a) Polarization of a GNR in a longitudinal electric field. Color scale of atom 
represents the intensity of induced charges. Dark arrows show the electric field around the 
graphcne. (b) Longitudinal polarizability per lnm width a*^ vs L. Circles denote simula- 
tion data and the dashed line stands for a numerically fitting curve. Coupled dependence 
of oijj on L and width d is shown in the inset, in which the color scale represents the 
amplitude of a/ / . 
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92 applied electric field, respectively. As shown in Fig. [3] (b), our simulation 

7/ 



93 results suggest that a), is roughly proportional to L 2 . 



L 2 

a*,, = — (L > Qnm) (4) 
B 

94 where B = 15.01 is a constant for either armchair or zigzag GNRs, since 

95 no large difference has been found between a*j of these two types of graphene. 

96 Since Eq. jljshows a typical metallic behavior of graphene, we note that a* j of 

97 semi-conducting graphene (minimum lateral dimension < 6nm [26j) should 

98 hold a linear relationship with L. Putting the empirical fits to a// and K* 

99 into Eq. |2j we finally obtain the governing equation of the electrostatic 

100 alignment of GNRs as follows: 

„ E 2 L sin [2(0 O -0)1 

6 = 11 (L > 20nm) (5) 

101 where C = 4AB pa 21eV-nm -2 . We note that, since the geometry periodic 

102 condition was applied in the width direction in our calculations, we effectively 

103 simulated graphene of infinite width. Thus, the edge disorder effects on 

104 structural [33] and mechanical properties [21] of graphene were neglected. 



105 4. Results and discussions 

we We can see from Eq. [5] that, in small deformation region, 9 is roughly 

107 proportional to E 2 , sin6 or L. In Fig. |4j we plot data of the alignment 

los ratio 6/6q as a function of E from molecular simulations in which periodic 

109 geometry condition is applied in width direction. A quantitative agreement 

no was obtained between the simulation data and the analytical prediction using 

in Eq. |5j The main inaccuracy can be considered to be contributed from the 
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Figure 4: (a) Cross sections of a GNR (L sa 9.2nm) in different electric fields E (9q = 7r/4), 
(a video recording file is available on line), (b) Alignment ratio 9/9q versus the field 
strength E for GNRs of different lengths. The symbols represent numerical data and the 
lines stand for results predicted by analytical model. 
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112 geometry approximations and the possible slight change of polarizability due 

113 to the curvature of graphene. We can see that the S-shaped curve of 9/9 
H4 tends to be flat near the maximal value 1 when the GNRs are well aligned 
115 to the field direction. In such a case that sin [2(9 — 9)] « 2(9 — 9), Eq. [5] 
lie can be simplified to 

9~ ~ 1 ~ C + 2E 2 L ^ 

H7 for large deformation. Compared with GNRs in same sizes, CNTs were 

us found to be less flexible in an electric field, e.g., a (5,5) SWCNT [L ~ lOnm) 

n9 can be bent to 9/9q ~ 0.4 in an electric field (E »s 2.0V/nm and 9q — 

120 while required field strength for producing the same amplitude of alignment 

121 for a GNR with the same width and length is about 7 times smaller (E w 

122 0.3V/nm). Analysis on the amplitudes of induced polarization in a GNR and 

123 a CNT shows that this difference is mainly due to the fact that the bending 

124 stiffness of the GNR is much lower than that of the CNT. Furthermore, 

125 we can predict that a multi-layered graphene should be less sensitive to an 

126 electric field than a single-layered one is, because of the electric screening 

127 effects and the friction between the layers [20J. 

128 5. Conclusion 

129 In conclusion, we have developed an analytical model to predict the align- 
no ment of cantilevered GNRs to an electrostatic field. Parameters used in this 

131 model such as the polarizability and the bending stiffness were determined 

132 from numerical fits to the data of simulations using an atomic electrostatic 

133 charge-dipole model and an empirical pseudo-chemical potential. This model 
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134 showed that the alignment angle roughly follows a linear relationship with 

135 the square of field strength and the graphene length when the deformation 

136 remains small. It was also found that, for GNRs with either armchair or 

137 zigzag edges, their effective bending stiffness and longitudinal polarizability 

138 are both approximately proportional to the square of graphene length. Com- 

139 parison showed that a GNR can be more easily oriented to electric fields than 
wo a CNT does, due to the GNR's low transverse bending stiffness. 
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